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Abstract 

Purpose Topsoil erosion due to land use has been 
characterised as one of the most damaging problems from 
the perspective of soil-resource depletion, changes in soil 
fertility and net soil productivity and damage to aquatic eco¬ 
systems. On-site environmental damage to topsoil by water 
erosion has begun to be considered in Life Cycle Assessment 
(LCA) within the context of ecosystem services. However, a 
framework for modelling soil erosion by water, addressing 
off-site deposition in surface water systems, to support life 
cycle inventory (LCI) modelling is still lacking. The objec¬ 
tives of this paper are to conduct an overview of existing 
methods addressing topsoil erosion issues in LCA and to 
develop a framework to support LCI modelling of topsoil 
erosion, transport and deposition in surface water systems, to 
establish a procedure for assessing the environmental damage 
from topsoil erosion on water ecosystems. 

Methods The main features of existing methods addressing 
topsoil erosion issues in LCA are analysed, particularly with 
respect to LCI and Life Cycle Impact Assessment methodol¬ 
ogies. An overview of nine topsoil erosion models is per¬ 
formed to estimate topsoil erosion by water, soil particle 
transport through the landscape and its in-stream deposition. 
The type of erosion evaluated by each of the models, as well 
as their applicable spatial scale, level of input data 
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requirements and operational complexity issues are consid¬ 
ered. The WATEM-SEDEM model is proposed as the most 
adequate to perform LCI erosion analysis. 

Results and discussion The definition of land use type, the 
area of assessment, spatial location and system boundaries are 
the main elements discussed. Depending on the defined sys¬ 
tem boundaries and the inherent routing network of the de¬ 
tached soil particles to the water systems, the solving of the 
multifunctionality of the system assumes particular relevance. 
Simplifications related to the spatial variability of the input 
data parameters are recommended. Finally, a sensitivity anal¬ 
ysis is recommended to evaluate the effects of the transport 
capacity coefficient in the LCI results. 

Conclusions The published LCA methods focus only on the 
changes of soil properties due to topsoil erosion by water. This 
study provides a simplified framework to perform an LCI of 
topsoil erosion by considering off-site deposition of eroded 
particles in surface water systems. The widespread use of the 
proposed framework would require the development of LCI 
erosion databases. The issues of topsoil erosion impact on 
aquatic biodiversity, including the development of character¬ 
isation factors, are now the subject of on-going research. 

Keywords Life cycle inventory • Surface water systems • 
Topsoil erosion 

1 Introduction 

Topsoil provides ecosystem services essential for life. It acts 
as a water depuration filter and as a substrate for growing 
food, fibre and biomass; it provides habitats for multiple 
organisms, contributing to biodiversity, and it provides a 
foundation for construction activities. Topsoil erosion is a 
natural and complex process that varies around the world 
depending on climate, land use, soil texture, ground slope, 
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vegetation cover, rainfall patterns, and field-level conservation 
practices (Montgomery 2007). Awareness of the on-site im¬ 
pacts of accelerated topsoil erosion began as early as the 1920s 
(e.g., Bennett and Chapline 1928), and since then, topsoil 
erosion has become a worldwide major environmental issue. 
This process occurs in three stages: detachment and entrain¬ 
ment of soil particles, their transport and eventual deposition. 
Topsoil erosion has impacts that are both on-site (the place 
where the topsoil is detached) and off-site (wherever the 
eroded soil particles reach the surface water systems). The 
eroded soil particles may come from several sources: soil 
erosion (by water, wind and tillage), mining and construction 
activities (Grimm et al. 2002). 

Because of the very slow rate of soil formation, any soil 
loss of more than 1 t ha -1 year -1 can be considered as irre¬ 
versible within a time span of 50-100 years. On average, the 
soil erosion rate for cropland is of the order of 7.51 ha -1 year -1 
(Grimm et al. 2002; Wilkinson 2005). Water is one of the 
major causes of soil erosion (Jones et al. 2012a), and the 
average soil erosion by water in Europe is about 
2.8 t ha -1 year -1 (Jones et al. 2012a, b) affecting 60 % of the 
total land area excluding Russia. However, losses as high as 
20^10 t ha -1 year -1 have been measured in individual storms 
in Europe (Jones et al. 2012a). 

In the last decades, the environmental damage to topsoil by 
water erosion has begun to be considered in Life Cycle As¬ 
sessment (LCA) within the context of ecosystem services (EC 
JRC 2010; Guinee et al. 2006). Special attention has been 
given to changes in soil properties at the local level, which 
may include loss of soil nutrients, salinisation, changes in soil 
organic matter, reduction of soil depth, increase of stoniness 
and reduction of water-holding capacity. 

However, there are also off-site impacts that affect aquatic 
biota caused by the eroded soil particles that are transported 
and deposited into surface water systems and into downstream 
soil. To date, these off-site impacts have not been considered 
in LCA, which is a substantial failing considering the scale of 
the problem. 

The existence of life cycle inventory (LCI) data on soil 
erosion by water is cmcial for the consideration of the associ¬ 
ated impacts in LCA. Currently, there is no regionalised and/ 
or local database of topsoil erosion by water. In addition, the 
spatial variability of soil texture and topographic parameters at 
a regionalised/local scale hampers the establishment of a 
conventional LCI. Therefore, to perform an unambiguous 
LCI of topsoil erosion, considering these spatial variabilities, 
it is necessary to use models based on geographic information 
systems (GIS). 

Several different models for predicting the topsoil erosion 
by water and the transport of the detached soil particles into 
the surface water systems, for different scales of catchment 
areas, have been developed (Kirkby et al. 2004; Merrit et al. 
2003; Van Rompaey et al. 2001a). However, these models 


have not been developed specifically for LCI purposes; they 
are applied primarily to watershed and river basin manage¬ 
ment issues, as well as to evaluate changes in water quality 
resulting from flood events and to assess siltation issues in 
catchment areas and navigable channels. 

In this study, in order to understand how topsoil erosion 
issues have been considered in LCA, an overview of existing 
methods addressing topsoil erosion is conducted. Further- 
more, to support LCI modelling and to fill the existing gap 
in current soil erosion research within the context of LCA, a 
framework for modelling soil erosion by water and deposition 
in water surface systems is proposed, firstly, an overview of 
the range of different models available for estimating the 
quantity of soil particles detached and entrained by water, 
their transport through the landscape and ultimate deposition 
into the surface water systems is performed. Subsequently, the 
topsoil erosion model that best fits the inventory of an LCA 
study is suggested. The constraints of the suggested model are 
underlined, and with the goal of operationalising the LCI 
modelling, some proper simplifications are recommended. 


2 Overview of LCA methods addressing topsoil erosion 

Because of the scale of soil erosion, the integration of topsoil 
erosion issues in the LCA structure assumes significant rele¬ 
vance, contributing to the establishment of environmental 
strategies for soil protection (e.g., CEC 2006). 

figure 1 illustrates the general cause-effect chain of topsoil 
erosion and shows the damage to human health and terrestrial 
and aquatic ecosystems. The topsoil erosion process causes 
local changes to soil properties and leads to off-site deposition 
of eroded soil particles in surface water systems and down¬ 
stream soil. This last pathway can bring off-site benefits for 
the downstream soil. When nutrients from the topsoil are 
transported from one place to another, the nutrients are 
redistributed downstream, which can lead to the increase of 
fertilisation of a downstream nutrient-deficient soil. 

Table 1 presents the methods developed for assessing the 
impact of topsoil erosion within the LCA context. All the 
methods developed by Cowell and Clift (2000), Mattsson 
et al. (2000), Muys and Garcia Quijano (2002), Garrigues 
et al. (2012) and Nunez et al. (2010) focus on topsoil erosion 
by water. However, Muys and Garcia Quijano (2002) ac¬ 
knowledge that when relevant, both wind and tillage erosion 
should be considered. For wind erosion, the authors suggest 
using wind erosion equations following the methodology of 
Schwab et al. (1993), whereas for tillage erosion, no recom¬ 
mendations are forthcoming. With regard to the level of 
damage, all methods analyse the local effects of topsoil ero¬ 
sion, i.e., the changes to the soil properties pathway (Fig. 1), 
focusing mainly on soil quality and soil depth indicators. 
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Fig-1 Cause-effect chain of topsoil erosion impacts. The shaded boxes relate to the pathway focused on in this study 


Nunez et al. (2013) applied the model following an end¬ 
point approach, whereas the models proposed by the other 
authors evaluate the effects of topsoil erosion at the midpoint 
level. At the midpoint level, these methods proposed different 
impact categories such as depletion of the soil resource, soil 
fertility and desertification depending on the method, whereas 
at the endpoint level, the area of protection from damage to the 
terrestrial ecosystem is assessed. Several of these methods are 
not dedicated exclusively to topsoil erosion damage. In Muys 
and Garcia Quijano (2002) and Garrigues et al. (2012), the 
indicator of the loss of soil is a sub-indicator in the soil sub¬ 
impact category; however, in Mattsson et al. (2000), the loss 
of soil is a sub-indicator in the soil fertility sub-impact cate¬ 
gory. These methods evaluate the topsoil erosion issues as a 
sub-category of the broader impact category of land use. In 
addition, at the endpoint level, Nunez et al. (2010) consider 
soil erosion as one of the four indicators (aridity, aquifer 
overexploitation, fire risk and erosion) for the desertification 
impact category. 

At the LCI phase, in the Cowell and Clift (2000), Muys and 
Garcia Quijano (2002) and Nunez et al. (2010) methods, the 
amount of topsoil eroded by water has been estimated by the 
Universal Soil Loss Equation (USLE) (Wischmeier and Smith 
1978). Garrigues et al. (2012) use the Revised Universal Soil 
Loss Erosion (RUSLE) equation (Renard et al. 1997). Both 
equations quantify the locally eroded soil particles but do not 
consider the runoff of these particles. However, the RUSLE 
equation presents some improvements in relation to the 
USLE, such as new calculation procedures to account for 
slope length and steepness, additional sub-factors for evaluat¬ 
ing the cover and management factors for cropland and 


rangeland and new conservation practice values for cropland 
and rangeland (Jones et al. 1996). Further details about the 
RUSLE equation are given in Section 4. 

With regard to the characterisation model, Cowell and Clift 
(2000) considered the Soil Static Reserve Life that is a func¬ 
tion of the global reserves of agricultural soil (total topsoil in 
the world, in t) and current annual global net loss of topsoil 
mass by erosion (in t year -1 ). Muys and Garcia Quijano 
(2002) recommend transforming the loss of soil mass into a 
loss of soil depth (in m) using the bulk density of the soil. 
Therefore, the loss of soil depth over a period of 100 years is 
compared with the total rootable soil depth up to 1 m. How¬ 
ever, neither method provides a set of operational characteri¬ 
sation factors (CFs). In fact, to derive such site-specific CFs, 
data on the reserve of the topsoil (area and depth of soil) under 
analysis is required. Both Mattsson et al. (2000) and Garrigues 
et al. (2012) use the loss of topsoil as an indicator for erosion 
impact. These authors argue that the estimates of loss of 
topsoil in the inventory level are informative enough to serve 
as impact indicators without requiring CFs. The characterisa¬ 
tion model suggested by Nunez et al. (2010) takes into ac¬ 
count the spatial area subjected to desertification, the decimal 
logarithm of the ecoregion area where desertification occurs 
and the time of perturbation. Moreover, this method proposes 
CFs for eight large natural areas, i.e., ecoregions (marine, 
prairie, temperate steppe, temperate desert, savanna, Mediter¬ 
ranean, tropical/sub tropical steppe and tropical/subtropical 
desert) that consider the loss of topsoil and terrain definition/ 
mass movement in each ecoregion soil weighted by the sur¬ 
face area of the ecoregion. Nunez et al. (2013) proposed an 
alternative characterisation model that considers the ultimate 
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Table 1 Overview of methods for inclusion of topsoil erosion issues in the LCA structure 
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damage to the depletion of the soil resource by assessing the 
loss of soil as the locally available soil reserves normalised 
with a reference soil depth of 3 m with transformation to 
emergy units (megajoule solar equivalents) per unit of area 
and time of land use. A set of regional CFs was developed at a 
spatial resolution of 10x10 km , taking into account the 
locally available soil reserves, the reference soil depth and 
the solar energy factor (the quantity of solar energy that is 
required to generate a gram of eroded soil particles) (Nunez 
et al. 2013). 

To our knowledge, there are no available LCI/Life Cycle 
Impact Assessment (LCIA) methods to assess the off-site 
effects due to topsoil erosion by water, i.e., the transport and 
deposition of eroded soil particles into surface water systems 
and the transport of soil particles from one place to another 
downstream (Fig. 1). Therefore, this study focuses on how to 
estimate topsoil erosion and to understand the transport of soil 
particles at the inventory level, evaluating the net quantity of 
soil particles that reaches the surface water systems. 

3 Models of topsoil erosion 

There are several purposes behind performing LCI modelling 
of topsoil erosion: 

• To describe as accurately as possible the topographic 
conditions of the land use system under analysis and the 
types of erosion that can occur (sheet, rill and gully 
erosion) 

• To characterise the soil profile, identifying the fraction of 
organic matter and the range size of the detached soil 
particles that reach the surface water systems 

• To understand the influence of tillage operations on the 
quantity of eroded soil 

• To quantify the quantity of detached soil particles that are 
transported and deposited into the surface water systems at 
sub-watershed level 

The topsoil erosion models that have become used most 
widely for predicting eroded soil particles and their transpor¬ 
tation and deposition into surface water systems are: 

• ANSWERS—Areal Nonpoint Source Watershed Envi¬ 
ronment Response Simulations (Beasley et al. 1980; 
Dillaha et al. 2001) 

• CREAMS—Chemical Runoff and Erosion from Agricul¬ 
tural Management System (Knisel 1980) 

• GUEST—Griffith University Erosion System Template 
(Yu et al. 1997; Rose et al. 1997) 

• HSPF—Hydrological Simulation Program Fortran 
(Johanson et al. 1980) 

• TOPOG (Gutteridge and Davey 1991) 


• WEPP—Watershed Erosion Prediction Project (Laflen 
et al. 1991) 

• MIKE-11 (Hanley et al. 1998) 

• WATEM-SEDEM—Water and Tillage Erosion Model- 
Soil particle Delivery model (Van Oost et al. 2000; Van 
Rompaey et al. 2001a) 

• EUROSEM—European Soil Erosion Model (Borselli and 
Torri 2010; Morgan et al. 1998) 

The goals, characteristics and limitations of each model are 
identified in Table 2. These models differ significantly regard¬ 
ing the type of erosion, their temporal and spatial scales, level of 
input data requirements, ease of use and in the algorithms used 
to predict the amount of detached and transported soil particles. 

The spatially distributed ANSWERS model was developed 
to evaluate the effect of agricultural land use practices on 
water quality of sub-watersheds. This model predicts the 
topsoil erosion by adopting the physics-based continuity 
equations used by Foster and Meyer (1972). The transport of 
soil particles and nutrients is estimated using a form of Yalin’s 
(1963) equation for bed load transport. The application of this 
model is limited by the large input data requirements, and it 
only considers sheet and rill erosion, disregarding gully ero¬ 
sion, which can be as extensive as the former two. Sheet 
erosion occurs by the detachment and entrainment of soil 
particles, which are then transported downslope by overland 
runoff (Hairsine and Rose 1992), whereas rill erosion occurs 
when the field slope increases, such that the overland runoff 
gains velocity, forming superficial channels (Rose 1993). 

The CREAMS model estimates the topsoil erosion and 
deposition of soil particles applying the USLE (Wischmeier 
and Smith 1978). The transport of soil particles is simulated in 
the overland flow using a steady-state continuity equation 
(Foster’s equation) at a plot sub-watershed level (range of 
4CM-00 ha). The unlikely assumption that all areas of the 
study are uniform in soil topography is a pertinent source of 
uncertainty and inaccuracy in the predicted topsoil erosion. 

The GUEST model predicts erosion, transport and deposi¬ 
tion using physics-based equations describing the steady-state 
soil particle flux (Rose 1993). The model was developed to be 
applicable at plot sub-watershed level (<100 ha). Despite 
being an adequate spatial scale to perform an LCI around 
agricultural or forestry land use, this model assumes uniform 
soil topography and does not consider tillage operations. 

In addition to topsoil erosion, both the HSPF and MIKE-11 
models allow the assessment of in-stream water quality pa¬ 
rameters (e.g., nitrogen, phosphorus, inorganic pollutants) at 
the watershed level. The HSPF model is applicable at the 
watershed level, in which the spatial area is divided into sub- 
regions with homogeneous hydrologic characteristics, i.e., 
homogenous edaphoclimatic data. The MIKE-11 model is 
also applicable at the watershed level, and the water system 
flow is described using physics-based St Venant equations. 
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Table 2 Goal, characteristics and limitations of models to predict topsoil erosion and the transport of soil particles to surface-water systems 


Models Goal and characteristics 


Limitations 


ANSWERS 


CREAMS 


GUEST 


HSPF 


TOPOG 


WEPP 


MIKE-11 


• Developed to evaluate the effects of agricultural land use practices 

in the on-stream water quality 

• Applicable at sub-watershed level (>400 ha) 

• A long-term oriented model using a continuous time step (e.g., 

days) 

• Considers sheet and rill erosion 

• Evaluates the effects of agricultural practices on pollutants in 

surface rainfall-runoff and topsoil erosion 

• Applicable at plot sub-watershed level (range of 40-400 ha) 

• Operates either on an event basis or long-term using a continuous 

time step (e.g., days) 

• Includes sheet, rill and ephemeral gully erosion 

• In addition to the prediction of topsoil erosion, also predicts the 

evapotranspiration of the crop, soil infiltration, among others 

• Developed to understand temporal fluctuations in soil particle 

concentration at the transport limit 

• Applicable at plot sub-watershed level (<100 ha) 

• An event-oriented model 

• Considers sheet and rill erosion 


• Developed to simulate the hydrology of watersheds and its in- 

stream water quality 

• Applicable at watershed level (~ 1,000 ha) 

• Operates either on an event basis or long-term using a continuous 

time step (from 1 min to 1 day, as long as the time step divides 
equally into 1 day) 

• Includes sheet, rill and ephemeral gully erosion 

• Takes into account tillage operations 

• Simulates the soil particle balance through the landscape until the 

reaching the surface water system 

• Simulates the transient hydrologic behaviour of watersheds and 

how this is affected by changing watershed vegetation and by the 
growth of the vegetation 

• Applicable at plot sub-watershed and watershed levels (range of 

10-1,000 ha) 

• Temporal resolution can be applied at daily time steps or sub-daily 

time steps 

• Uses detailed topographic information that can be provided by a 

digital elevation model of the area of study 

• Considers sheet and rill erosion 

• Developed to estimate net soil loss, the effect of hydraulic 

structures and impoundments on runoff flow and soil particle 
transport. 

• There are two versions: one that is applicable at plot sub-watershed 

(< 100 ha) and a second applicable at watershed level 
(-1,000 ha) 

• Operates either on an event basis or long-term using a continuous 

daily step 

• Includes sheet, rill and ephemeral gully erosion 

• Takes into account tillage operations 

• Used to predict topsoil erosion and to model the in-stream water 

quality 

• Applicable at watershed level 

• Operates either on an event or long-term basis 

• Uses detailed topographic information that can be provided by a 

digital elevation model of the area of study 

• Considers sheet and rill erosion 


• Disregards gully erosion 

• Complex model requiring large input data (e.g., river network, total 

porosity of soil, infiltration control zone depth, erodibility, steady 
state infiltration, among others) 

• Does not take into account tillage operations 

• Requires calibration and validation 

• Predicts topsoil erosion following the one-dimensional USLE 

approach 

• Assumes uniform topography (slope) of soil and land use 

• Does not take into account tillage operations 

• Complex model with large input data requirements (e.g., 

precipitation series, monthly air temperature and solar radiation 
values, crop type data, among others) 

• Assumes uniform topography (slope) of soil 

• Veiy complex model requiring large input data (e.g., runoff rate, 

length, width and slope of the soil, percent of sand grains of 
primary particles, soil particle/water-stable aggregate size 
distribution, among others) 

• Does not take into account tillage operations 

• The general lack of input data means that considerable effort is 

required to use it to predict topsoil erosion without prior 
calibration against parameterisation of field measurements 

• As inputs requires of precipitation data, estimates of potential 

evapotranspiration, topography, solar radiation, humidity, among 
others 

• Relies greatly on calibration against parameterisation of field 

measurements 


• Disregards gully erosion 

• Does not take into account tillage operations 

• Requires calibration 


• Large computational and input data requirements (e.g., leaf area 

index, canopy cover and height, duration of runoff, bulk density 
of the soil, among others) 

• Relies on calibration against parameterisation of the field 

measurements 


• Large input data requirements (e.g., hydrometric and topographic 

data, moisture content in surface and root zone, among others) 

• Does not take into account tillage operations 

• Low-quality data of the input parameters 

• Relies greatly on calibration against parameterisation of the field 

measurements 
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Table 2 (continued) 


Models 

Goal and characteristics 

Limitations 

WATEM- 

• Used to simulate eroded soil particles by the RUSLE method 

• Low quality data of the input parameters 


SEDEM 


EUROSEM 


considering a two-dimensional landscape structure and the 
transport and deposition of soil particles into surface water 
systems 

Applicable at plot sub-watershed and watershed levels 
A long-term average annual deposition of soil particles in the 
surface-water systems 

Focuses on the spatial resolution of 20x20 m 
Includes sheet, rill and ephemeral gully erosion 
Uses detailed topographic information that can be provided by a 
digital elevation model of the area of study 
Takes into account tillage operations 
Requires a modest amount of input data parameters 
The parameter-based nature of the model allows an easy analysis of 
the contribution of individual parameters (sensitivity analysis) to 
the predicted topsoil erosion and deposition 
Predicts topsoil erosion 

Applicable at plot sub-watershed level (<100 ha) 

Operates on an event basis considering a temporal resolution of 
minute-by-minute 


Requires partial calibration on the soil particle transport capacity 


Only considers sheet and rill erosion 
Large input data requirements (rainfall data, soil mechanical 
properties, topographical information, micro-topographical 
information, soil properties and vegetation information) 
Does not take into account tillage operations 
Low-quality data of the input parameters 
Requires calibration and validation 


Both models require large amounts of input data and over¬ 
look the relevance of gully erosion in a similar way to the 
ANSWERS and GUEST models. Gully erosion describes 
deep channels of runoff water and soil particles. 

The TOPOG model describes water movement through the 
landscape, over the land surface, into the soil, through the soil 
and groundwater and back to the atmosphere via evaporation. 
Topsoil erosion is simulated, and the transport of soil particles 
is simulated following the Engelhund and Hansen (1968) 
equation. It disregards gully erosion processes, and it requires 
large and detailed input data. 

The WEPP model evaluates net soil loss as well as the 
effects of hydraulic structures and impoundments on runoff 
flow and soil particle transport. It assesses the effects of tillage 
operations on the quantity of soil particles deposited into the 
surface water systems. The erosion and transport components 
are determined using the same relationships as the AN¬ 
SWERS model. It tends to overpredict the transport of soil 
particles for small events and underpredict their transport for 
large events (Yu 2005). Gully erosion is disregarded at plot 
sub-watershed level, which does not happen at the watershed 
level. 

The WATEM-SEDEM model predicts the long-term mean 
annual topsoil erosion following the RUSLE equation 
(Renard et al. 1997). It is applicable at plot sub-watershed 
and watershed levels, and it allows the evaluation of the 
influence of tillage operations on the quantity of eroded soil. 
Sheet, rill and ephemeral gully erosion processes are 
considered. 


The EUROSEM model predicts topsoil erosion and simu¬ 
lates soil particle transport to the surface water systems at the 
plot sub-watershed level. However, it only simulates single 
intense rainfall-runoff events, i.e., long-term simulation of 
multiple rainfall-runoff events is a performance constraint. It 
overlooks the relevance of gully erosion in a similar way to the 
ANSWERS, GUEST, HSPF, MIKE-11 and TOPOG models. 

4 Life cycle inventory modelling of topsoil erosion 
and deposition 

4.1 Model suggestion and its overall structure 

Because of the local or regional nature of topsoil erosion, the 
spatial differentiation assumes particular relevance. A model 
that performs an erosion-deposition inventory at a plot scale 
and uses site-dependent parameters allows the reduction of the 
uncertainty of the erosion LCI compared with generic default 
data on a broader scale. The ANSWERS, HSPF and MIKE-11 
models have been developed to be applicable on a broad 
spatial scale. In addition, these models present operational 
complexity, require large amounts of input data, evaluate in- 
stream water quality, which is not an aim of the soil-deposition 
inventory, and require calibration. Although the ANSWERS 
model considers gully erosion and allows the inclusion of 
tillage operations, it is only applicable at the watershed level, 
i.e., the scale of assessment is approximately 1,000 ha, hin¬ 
dering the topsoil analysis of a local land use system. For these 
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reasons, these models are not adequate to perform an LCI of 
an agricultural or forestry system on a plot scale. The 
CREAMS and GUEST models also do not fit the erosion- 
deposition LCI purposes. Despite being applicable at a local¬ 
ised scale, these models consider uniform soil topography. 
Any model that uses a digital elevation model (DEM) to 
define the two-dimensional slope length parameter will reduce 
the underestimation of the erosion-deposition inventory of a 
non-uniform soil and land use. The shape of a slope affects the 
average topsoil eroded and the transport along the hill slope to 
the surface water system. For instance, according to 
Wischmeier and Smith (1978), the eroded soil from a convex 
slope can easily be 30 % greater than that from a uniform 
slope. In addition, these models disregard gully erosion and 
need calibration. 

The EUROSEM model also does not fit the LCI purpose. It 
operates on an event basis, i.e., by each rainfall-runoff event. 
This means that to estimate the average eroded soil particles 
transported to the surface water systems during a sown- 
harvested crop land use study or forestry land use study, 
numerous simulations should be performed to obtain the 
cumulative erosion-deposition inventory. In addition, gully 
erosion is not considered, and it requires a large amount of 
input data, as mentioned in the previous section. 

The WEPP model could be adequate for erosion-deposition 
inventory purposes as it is applicable at a plot sub-watershed 
level, operates on a long-term basis, considers gully erosion 
and takes into account the influence of tillage operations in the 
erosion and transportation of soil particles. However, the large 
computational and input data requirements and the calibration 
needed for many of the model parameters (e.g., slope, land 
cover, soil type) are shortcomings preventing its widespread 
use in LCI studies. 


The WATEM-SEDEM model has the same characteristics 
as the WEPP model but requires fewer input data parameters. 
It is the only model that uses RUSLE (Renard et al. 1997) in 
considering the variability of soil texture and topographic 
parameters (elevation, slope, upslope, profile curvature). This 
acquires cmcial relevance in establishing the routing network 
transport of soil particles to the water systems. Indeed, the 
primitive form of this equation, the USLE (Wischmeier and 
Smith 1978), and its update, the RUSLE, have been well- 
accepted by the LCA community for estimating the annual 
mass of topsoil loss, as mentioned in Section 2. However, both 
the primitive and RUSLE equations account only for the 
detached and entrained stages of the soil particles. The signif¬ 
icant variability of soil texture and topographic parameters 
suggests the need for GIS to account for a less crude and 
unambiguous inventory of the detached soil particles deposit¬ 
ed into the surface water systems. Indeed, performing a 
geospatial LCI is a step forward from conventional LCI stud¬ 
ies. After considering the main characteristics and constraints 
of each model regarding performing an LCI of topsoil erosion, 
presented in Section 2, it is suggested that the use of the 
WATEM-SEDEM model is most appropriate. 

The WATEM-SEDEM is a simple topography-driven model 
(Fig. 2) that uses the RUSLE equation (Eq. (1)) to estimate the 
topsoil eroded by water per land area (Renard et al. 1997): 

A = R x K x LS 2D xC xP (1) 

_o _i 

where A is the average eroded soil per year (kg m year ), R 
is the rainfall-runoff erosivity parameter 
(MJ mm m h year ), K is the soil erodibility parameter 
(kg h MU 1 mm -1 ) that represents the soil resistance to erosion 
by water and is related to the soil characteristics, namely, the 


Inputs - topsoil erosion 


Land use map: 

identifying the land use 
categories: cultivated 

land, forest, grassland, 
rural areas and water 
systems 


Digital elevation model 

(DEM) to calculate the 
LS2D 


Crop management 
factor (C) for each land 
use category: cultivated 
land, forest and grassland 


Erosivity parameter (R) 

calculated from monthly 
rainfall data 

- 

Erodibility parameter 

(K): calculated from 
organic matter, texture, 
structure, permeability, 
or collected as secondary 
datum 


r 



Inputs - tillage operations 


r > 

Tillage transport 
coefficient (k till ) 




J 




Eroded topsoil 


Inputs - transport of soil particles 


The same inputs + 


Transport capacity 
coefficient (ktc) for 
non-erodible land 
surfaces (ktc min ) and for 
arable land (ktc max ) 



Soil particles transported 
and deposited in the 
surface water systems 


7 1 

Fig. 2 An overview of the WATEM-SEDEM model structure distinguishing the topsoil eroded (kg m year ), the transport of soil particles by 
overland runoff and their deposition in surface water systems (kg m 1 yeah 1 ) 
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structure, texture, organic matter content and permeability, 
LS 2 d is the two-dimensional slope-length parameter 
(dimensionless), C is the cover-management parameter 
(dimensionless) allowing an understanding of the extent to 
which vegetation cover prevents soil erosion, and P is the 
support-practice parameter to reduce runoff and soil erosion 
(dimensionless). 

The tillage erosion process is also simulated by the 
WATEM-SEDEM model following a diffusion-type equation 
(Eq. (2)) developed by Govers et al. (1994): 

= (2) 

where Q st is the rate of net downslope soil transport per tillage 
translocation (kg m 1 per tillage operation), k ti \\ is the tillage 
transport coefficient (kg m -1 per tillage operation), h is the 
height at a given point of the hill slope (m), and x is the 
distance in the horizontal direction (m). 

Unlike water erosion, topsoil displacement will only occur 
within the field, i.e., it is not transported directly to the surface 
water systems. However, these detached soil particles can 
contribute to the increase in the quantity of soil particles that 
are transported by overland runoff flow. This means that soil 
tillage operations have a non-negligible influence on topsoil 
erosion (Govers et al. 1994), and thus, they should be included 
in an LCI study. 

The transport of the detached soil by overland runoff to the 
surface water systems is calculated according to Eq. (3). It is 
assumed that the transport capacity is proportional to the 
potential gully erosion (Desmet and Govers 1995; Van Oost 
et al. 2000; Van Rompaey et al. 2001a): 

Tc = ktcxtfx^x (LS 2D -4.12 x Sg 0 ' 8 ) (3) 

where Tc is the transport capacity (kg m -1 year -1 ), ktc is the 
transport capacity coefficient (m), R , K and LS 2 d are the 
parameters from the RUSLE equation, and Sg is the local 
slope gradient (dimensionless). 

4.2 Guidelines and simplifications 

In this section, the guidelines, methodological simplifications 
and recommendations for the LCI modelling of topsoil ero¬ 
sion using the WATEM-SEDEM model are identified and 
described. 


4.2.1 Geographical area and multifunctionality 

To estimate the quantity of detached soil particles from a land 
use system that reach the surface water systems, the geograph¬ 
ical land area should be defined. Moreover, the LCI should 
include the classification of the land use type under analysis. 


After the definition of the geographical land area, the next 
step is to delineate the related DEM, which represents the 
topographical variability of the soil. In addition to the GIS 
information, it is also necessary to establish a soil particle 
routing network, which represents the runoff flow through 
the landscape towards the surface water systems. For this, a 
land use map identifying the different land use types (for 
instance, cultivated land, forest, grassland, rural areas and 
water system) is required as an input to the WATEM- 
SEDEM model (Alatorre et al. 2010). The establishment of 
the soil particle routing network can raise the problem of 
incoherence between the pre-defmed geographical land area 
under study and the overall downstream landscape, because 
the delineated DEM of the land use system under study 
may not match the land use map. To overcome this 
constraint, the DEM should be expanded to include all 
the downstream land area towards the surface water 
system. This highlights a system multifunctionality 
problem, because an overestimated erosion-deposition 
inventory is performed, i.e., it includes soil particles 
detached and transported from the land use system 
under study and downstream land use systems. As a 
result, in most cases, it is not possible to conduct an 
LCI study through a single run of the WATEM-SEDEM 
model. In order to assign the quantity of detached soil 
particles from the pre-defmed land use system under 
study and its delivery to the surface water systems, a 
system subdivision should be carried out, as illustrated 
in Fig. 3. By considering the agricultural/forest land 
area under study and the downstream land area routing 
to the surface water system, two co-functions are pro¬ 
vided by this unique system (AB) (the topsoil erosion 
process from the system under study and the topsoil 



System AB 

Agricultural/forestry land area under study + downstream 
land area routing to the surface water system 


System A 

Agricultural/forestry land area 
under study 




Soil particles routing' Surface water system 



Soil particles from system B 


Soil particles from system A and B 




Fig. 3 Solving the system multifunctionality by subdivision of system 
AB. Subdivision aims to obtain the LCI topsoil erosion of system A 
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erosion process from the downstream landscape). Fig¬ 
ure 3 shows that the inventory of the pre-defmed land 
use system (system A) results from the subdivision of 
the overall system (AB). In practice, the LCI practition¬ 
er must have in mind that the spatially distributed 
WATEM-SEDEM model needs to be run at least twice 
to establish the inventory: the first run accounts for the 
inventory of the overall system (AB) and the second 
run accounts for the inventory of system B. 

The exception occurs when the pre-defmed land use 
system is located in the water system bank. In this case, 
it is reasonable to assume that all the eroded soil is 
transported to the water. In this case, the erosion- 
deposition inventory can be performed by the direct 
use of the RUSLE equation. 


4.2.2 Input data requirements and simplifications 

The input data requirements to perform an LCI model¬ 
ling using the WATEM-SEDEM model are the soil 
erodibility parameter, rainfall-runoff erosivity parameter, 
two-dimensional topographical parameter, cover manage¬ 
ment parameter, tillage transport coefficient and the 
transport capacity coefficient of the soil particles. 

Soil erodibility parameter (K) The soil erodibility parameter 
( K ), as mentioned previously, describes the susceptibil¬ 
ity to erosion by rainfall, and it depends on the texture, 
structure, organic matter, permeability and physico¬ 
chemical interactions of the soil. Generally, the avail¬ 
ability of detailed primary soil data and soil erodibility 
parameters at the local scale is scarce. Therefore, in¬ 
stead of a soil erodibility map with a 20x20-m spatial 
resolution, the LCI practitioner can assume that the 
overall system (AB) has a single and a long-term average 
K parameter. 

Depending on the soil profile, the average K param¬ 
eter can be estimated following different approaches. If 
the silt fraction of soil does not exceed 70 %, the K 
value (kg h MJ _1 mm" 1 ) can be estimated using the 
nomograph developed by Wischmeier and Smith (1978) 
or by applying its algebraic approximation (Eq. (4)) 
(Renard et al. 1997): 

2.1 x 10~ 4 (12-OM)M U4 +3.25(5-2) + 2.5(p-3) 

“ 100 x 7.59 

( 4 ) 

where OM is the fraction of organic matter (%), M is the 
particle-size parameter (percentage silt+fme sand fraction 
content) x (100-clay fraction), and s and p are the soil structure 
(dimensionless) and permeability classes (dimensionless), 
respectively. 


If the soil profile has more than 70 % silt, the K parameter can 
be calculated following Eq. (5) (Wischmeier and Smith 1978): 


K 


7.594< 0.0034 + 0.0405exp 


-0.5 


/logD g + 1.659 
V 0.7101 



where D g is the geometric mean weight diameter of the 
primary soil particles (mm), which can be measured experi¬ 
mentally by analysing soil samples. This relationship is only 
validated for soil with less than 10 % of rock fragments by 
weight (fraction higher than 2 mm). When this equation is not 
valid, it is recommended to collect the input K parameter 
datum from literature (e.g., Borselli et al. 2009; Gomez et al. 
2003; Irvem et al. 2007; Sanchis et al. 2008), national or 
international research projects (e.g., Van der Knijff et al. 
2000) and soil databases attending to the texture of the soil 
under study. For instance, depending on the soil profile, this 
secondary datum can be collected from the Land Use and 
Cover Area frame Survey database (Panagos et al. 2012). 


Rainfall-runoff erosivity parameter (R) The rainfall-runoff 
erosivity parameter ( R ) represents the impact of rain on topsoil 
erosion, and it is assessed based on monthly rainfall data that 
can be collected from meteorological stations within the land 
area under study and for a long-term period (for instance 20- 
30 years). When there are no meteorological stations within 
the land area, the average monthly precipitation data can be 
taken from the NewLocClim model (FAO 2005), which 
provides estimations of average climatic conditions at loca¬ 
tions for which no observations are available, considering the 
nearest neighbour interpolation method. In these situations, 
the R parameter can be calculated using Eqs. (6), (7) and (8) 
(Renard and Freimund 1994): 



R = 0.7397 x MF 1 ' 847 , for MF < 55mm (7) 


R = 95.77-6.081 x MF + 0.4770 x MF 2 , forMF > 55mm 

( 8 ) 

where MF is the Modified Fournier index (mm) (Amoldus 
1977), is the average monthly precipitation (mm), and P is 
the average annual precipitation (mm). Attending to the area 
of assessment, i.e., plot sub-watershed level, it is recommend¬ 
ed to assume that the spatial variability of parameter R is 
negligible. Therefore, the long-term variability of R on a 
yearly basis is considered uniform for the entire study area. 

Two-dimensional topographical parameter (LS 2 d) The two- 
dimensional topographical parameter (LS 2 d) shows the 
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impact of the length and slope of the landscape on soil erosion, 
and it has a large value whenever the length and slope of the 
landscape are large. In WATEM-SEDEM, the DEM is used to 
calculate the LS 2 d parameter. 

Cover management parameter (C) The cover management 
parameter (Q reflects the effect of cropping and management 
practices on erosion rates, i.e., it indicates the extent to which 
the vegetation cover prevents topsoil erosion. This parameter 
depends on the size of cover plants, the state of the surface 
area, plant roots, surface roughness and amount of contained 
water (Park et al. 2011). There is very limited availability of C 
parameters for different crops and/or forests at different loca¬ 
tions (Nunez et al. 2013). Because of this lack of data, Nunez 
et al. (2013) developed a methodology to estimate specific C 
parameters for agricultural and/or forestry systems. However, 
this approach requires a follow-up field to calculate the per¬ 
centage of vegetative cover throughout the sown-harvest pro¬ 
cess. To overcome this time-consuming procedure, we suggest 
using single values concerning the land use categories included 
in the defined overall system (AB), as well as in system B: crop, 
forestry and/or grassland. Even though this increases the un¬ 
certainty of the inventory results, it is the best way in which to 
minimise the complexity of calculating C parameters at the plot 
scale. If the system (AB) encompasses various types of crops 
and/or forestry land use, a weighted average between each 
culture cover management and the occupied area should be 
performed. The C parameters can be collected from scientific 
publications (e.g., Lee and Lee 2006; Keesstra et al. 2009; Yang 
et al. 2003), books and reports (e.g., Almorox et al. 1994; 
Pimenta 2003; Van der Knijff et al. 2000). 


2006) and in Spain (Alatorre et al. 2010); values ranging 
between 1:3.80 and 1:2.20 for mountainous and non- 
mountainous areas in Italy (Van Rompaey et al. 2005); 
1:3.89 for seven small catchments in South Africa (Van 
Rompaey et al. 2001b; Verstraeten et al. 2001) and 1:2.50 
for the Czech Republic (Van Rompaye et al. 2003). The LCI 
practitioner should compare the topographic characteristics 
(e.g., LS 2 d) of the land use system under study with the 
characteristics of the several regions mentioned above. Re¬ 
gions with similar topographic characteristics should be 
characterised by similar ktc max /ktc min ratios. 

Data on annual suspended soil particles input into the water 
system channels under study, as well as data on rainfall 
erosivity for the same time series, are required to perform 
the calibration. For each combination of ktc max and ktc min , a 
quantity of soil particles delivered into the surface water 
system is predicted, allowing the comparison of the quantities 
predicted by the model and those measured. The Nash- 
Sutcliffe model efficiency statistic (NS) (Nash and Sutcliffe 
1970) is used as a measure of likelihood following Eq. (9): 


E XOrPti 

NS = l- ^ i=1 - 

el, m: 



where n is the number of observations, Oi is the measured soil 
particles in the water system, P t is the predicted soil particles 
in the water system, and O is the average of the measured soil 
particles. NS can range from -oo to 1 and represents the 
proportion of the initial variance accounted for by the model. 
The closer the value of NS is to 1, the more efficient is the 
model (Alatorre et al. 2010). 


Transport capacity coefficient of soil particles (ktc) The trans¬ 
port capacity coefficient of soil particles (ktc) depends on the 
local land vegetative cover and represents the slope length 
needed to runoff an equivalent quantity of topsoil particles 
from a bare surface with a similar slope gradient (Verstraeten 
2006). The WATEM-SEDEM model requires distinction be¬ 
tween arable land surfaces, ktc max (cultivated land), and non- 
erodible land surfaces, ktc min (forests and/or grassland). More¬ 
over, Van Rompaey et al. (2001a) have indicated that a cali¬ 
bration of transport capacity coefficient would be desirable 
because of the unique and specific routing network of soil 
particles for each land use system. In fact, land use systems 
with high rainfall erosivity also present high runoff flows and 
soil particle deposition into the surface water system resulting 
in higher transport capacity (Verstraeten et al. 2007). Because 
of the unavailability of sufficient data to calibrate ktc max and 
ktc min separately, Verstraeten (2006) suggested maintaining a 
fixed ratio of ktc max /ktc min values. For instance, the following 
typical values between ktc min and ktc max have been used: 
1:3.33 in central Belgium and the southwestern part of Slove¬ 
nia (Keesstra et al. 2009; Verstraeten et al. 2006; Verstraeten 


Tillage transport coefficient (k ti u) The intensity of topsoil 
erosion is influenced by the tillage transport coefficient (£ t m)- 
Experimental tillage erosion studies have shown that tillage 
erosivity is affected by parameters other than slope gradient, 
such as tillage depth, tillage speed and soil condition 
(Gerontidis et al. 2001). Van Muysen et al. (2000) presented 
a dataset of (£ ti n) calculated for different implements, tillage 
speed and tillage depth. In addition, (& ti n) following the up- 
and down-slope tillage direction for each implementation can 
be calculated following Eq. (10) (Van Muysen et al. 2002): 

k m = 2.026 x p b x D 1 ’ 989 x v 0 406 (10) 

_O 

where p b is the bulk density of the soil (kg m ), D is the 
tillage depth (m), and v is the tillage speed (m s -1 ). For a set of 
tillage operations during the sown-harvest cycle, the tillage 
transport coefficient is the sum of the coefficients for each 
individual operation. 

When the system (AB) encompasses cultures that require 
different tillage operations, a weighted average tillage trans¬ 
port coefficient should be used. 
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4.2.3 Complexity ofLCI modelling and uncertainty 
parameters 

When defining the scope of the study, special attention should 
be given to the definition of the geographic system boundary, 
identifying and reporting the location of the land use system. 
On the one hand, depending on the spatial location of the land 
use, different modelling procedures can be followed. In some 
situations, as described in Section 4.2.1, a system subdivision 
is required to solve the multifunctional and to predict the 
erosion-deposition inventory of the pre-defmed land use 
system. 

On the other hand, depending on the area of assessment, it 
can include different land use types and huge variability of soil 
properties. From an ideal perspective, a map that could show 
the variability of each input parameter necessary to estimate 
the quantity of topsoil eroded would be desirable. 

However, currently, this is an unrealistic procedure because 
of the non-availability of the required input data. First, there 
are no specific or site-generic LCI databases containing soil 
erodibility and cover management maps. In particular, data 
related to cover management parameters for specific cultures 
and locations are quite rare. Secondly, although it is possible 
to perform experimental measurements during the overall 
growth of the crop/forestry from sowing to harvest, this is 
impracticable from a time and economic point of view. There¬ 
fore, to overcome these constraints and to operationalise the 
LCI inventory, it seems reasonable to use an average value for 
each input parameter, i.e., constant values instead of maps. 

The overall uncertainty of the erosion-deposition LCI is 
caused by the empirical inaccuracy of the input data, unrep¬ 
resentative data, lack of data sets and inherent uncertainties of 
the model. Perhaps the biggest problem with the modelling of 
eroded soil and its transport to the surface water systems, 
related to the availability and reliability of input data, is that 
the transport capacity coefficient is a major source of 
uncertainty and inaccuracy in an LCI performed with 
the WATEM-SEDEM model (Van Rompaey et al. 2001a; 
Verstraeten et al. 2006). To understand the influence of a set 
of different ktc max /ktc min ratios on the mass of soil particles 
detached and deposited, a sensitivity analysis is recommend¬ 
ed. Moreover, the influence of other input parameters, for 
instance, K , C and k t m , which are considered constant and/ 
or adapted from literature, should also be evaluated. 

5 Conclusions 

The existing LCA methods addressing soil erosion focus only 
on the local changes of soil properties. They differ in the way 
that they quantify erosion at the LCI level and in the LCIA 
methodology, namely, in the assessment level (midpoint/end- 
point), impact pathway, category indicator, characterisation 


model and characterisation factors. None of these methods 
considers the transport and deposition of the eroded particles 
towards the surface water systems. This study provides the 
first step for the inclusion of this issue in the LCA structure. 
The conducted overview of erosion models shows that there is 
a range of available models for predicting the quantity of 
detached soil particles and their transport towards the surface 
water systems. Following consideration of their operational 
complexity, input requirements and spatial scales, it has been 
suggested that the WATEM-SEDEM model is most appropri¬ 
ate for performing an LCI analysis. Therefore, this paper 
provides a framework with which to perform an LCI of the 
topsoil eroded and transported towards the surface water 
systems. This operational framework has the following key 
points: 

• Definition of the spatial land area of assessment (system 
boundary) stating its geographical location. 

• When necessary, a system subdivision should be per¬ 
formed to avoid allocation procedures, depending on the 
pre-defmed land use systems (allocation procedures). 

• Collecting all the required input data to run the WATEM- 
SEDEM model. 

• Performing a sensitivity analysis to evaluate the influence 
of the transport capacity coefficient on the LCI results. 

This framework seems to be a reasonable approach to 
establish realistic estimations of the order of magnitude of 
topsoil erosion and the amount of soil particles that reach the 
surface water systems, causing turbidity and in-stream soil 
particle deposition, which in turn affect aquatic biota. 

It should be noted that this framework should be 
operationalised by performing case studies. Moreover, 
displaced soil particles are themselves a source of potential 
environmental harm, especially when they reach water sys¬ 
tems. Therefore, a method to evaluate the impact of topsoil 
erosion on aquatic biodiversity within the LCIA is also re¬ 
quired, and it is the subject of on-going research. 
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